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One fascinating way of revealing quantum nonlocality is the all-versus-nothing test due to Green- 
berger, Home, and Zeilinger (GHZ) known as GHZ paradox. So far genuine multipartite and mul- 
tilevel GHZ paradoxes are known to exist only in systems containing an odd number of particles. 
Here we shall construct GHZ paradoxes for an arbitrary number (greater than 3) of particles with 
the help of qudit graph states on a special kind of graphs, called GHZ graphs. Furthermore, based 
on the GHZ paradox arising from a GHZ graph, we derive a Bell inequality with two d-outcome 
observables for each observer, whose maximal violation attained by the corresponding graph state, 
and a Kochen-Specker inequality testing the quantum contextuality in a state-independent fashion. 



Local realism cannot make quantum theory complete, 
as argued by Einstein, Podolsky, and Rosen (EPR) based 
on the belief that every element of physical reality must 
have a counterpart in a complete theory 1 . Accord- 
ing to them, an element of reality is corresponding to 
a physical quantity whose value can be predicted with 
certainty without in any way disturbing a system. No 
disturbance is ensured by the locality, i.e., the assump- 
tion that the result of a measurement cannot be affected 
by any spacelike separated events. The clashing between 
the local realism and quantum mechanics as revealed 
by several no-go theorems such as Bell's theorem [2J, 
Greenberger-Horne-Zeilinger (GHZ) theorem [MS], and 
Kochen-Specker (KS) theorem [6], shows that the quan- 
tum mechanical description of our world is nonlocal, or 
more generally contextual. This fascinating and funda- 
mental quantum feature of nonlocality and contextual- 
ity has been verified in experiments on various physical 
systems, e.g., [7], via the detection of violations of Bell 
inequalities and KS inequalities [SrtTU]. 

Among these genius approaches, GHZ theorem [3] H] 
provides us an "all-versus-nothing" [IT] test of a stronger 
type nonlocality, referred to as GHZ nonlocality, than 
Bell's nonlocality. This is a state-dependent argu- 
ment: because of the perfect correlations in some spe- 
cial state called the GHZ state, e.g., a 3-qubit GHZ 
state |$) = ^(|000) — |1H)), some local observables 
are elements of reality according to EPR. For example 
since the observable a^a^ay stabilizes the GHZ state, 
i.e., a],o-y(jy |$) = |$), observables cr;). , a 2 , cr^ are all el- 
ements of reality. Here cr^ denote 3 standard Pauli 
matrices for the fcth qubit. Similarly, from two other 
stabilizers a y a^ay and a^a^a^ of |$) we know that all 
a* are elements of reality and must have realistic val- 
ues m,u' v = ±1 for k = 1,2,3. Realistic values are sup- 
posed to obey the same algebraic relations as their cor- 
responding observables. That is to say we have on the 
one hand mfmfmf = —1, since a^cr 2 o^|<I>) = — 1$) and 
m^m|m| =mj , m|m| =m v 1 rn\m : ^ — 1. On the other 



hand, since {ra v k ) 2 = 1, we have identity mfmfmf = 
(mf m\m\){m\m2m\){m\my 2 rri2 ) which gives rise to a 
contradiction —1 = 1. 

This elegant presentation of the GHZ paradox for 3 
qubits is due to Mermin [3] soon after its first discovery 
for a 4-qubit GHZ state [3] and has already been ver- 
ified experimentally [12]. Although originally the GHZ 
argument is state dependent, it was found recently that 
any GHZ paradox can give rise to a KS inequality for 
a state- independent test of quantum contextuality [H]. 
In addition to its fundamental role played in our under- 
standing of quantum nonlocality and contexuality, the 
GHZ paradox also finds numerous applications such as 
in the quantum protocols for reducing communication 
complexity [13] and for secret sharing [13] . 

Compared to the bipartite and two-level case, mul- 
tipartite and multilevel nonlocality or entanglement is 
poorly understood. In some quantum informational tasks 
such as quantum cryptography, the usage of multidimen- 
sional systems offers advantages such as an increased level 
of tolerance to noise at a given level of security and a 
higher flux of information compared to the two dimen- 
sional case[J5]. Thus, it is crucial to investigate the rel- 
evant physical properties from some subclasses of these 
systems, e.g., GHZ nonlocality from a special kind of qu- 
dit states. Earlier efforts (TB][T7] to generalize GHZ para- 
doxes to multidimensional and multilevel systems can be 
reduced either to the qubit cases or to fewer particle 
cases, except the cases of n — 4j + 3 for qubits [17 1. Gen- 
uine multipartite multilevel GHZ paradoxes were first 
found by Cerf et al. for (d + l)-partite d-level systems 
with d being even [18] . An unconventional approach by 
using concurrent observables, not commuting yet hav- 
ing a common eigenstate, is proposed by Lee et al. to 
construct a GHZ paradox for the GHZ states of an odd 
number of particles [15] ■ Also a GHZ-like argument (all- 
versus-something) is proposed by Kaszlikowski et al. for 
c?-partite d- level systems [20] , in which concurrent observ- 
ables have been used implicitly. Later, DiVinccnzo and 
Peres [2T| found out that not only can GHZ states exhibit 
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the GHZ paradox but also those code words, which are 
one kind of multipartite entangled states used in quan- 
tum error corrections [22], can exhibit GHZ nonlocality. 
But so far genuine multipartite and multilevel GHZ para- 
doxes for an even number of particles are still missing. 

It turns out that GHZ states as well as code words 
from stabilizer codes [53] are graph states [23] which are 
essential resources for the one-way computing |25j and 
also provide an efficient construction of quantum error- 
correcting codes [2S]. It is thus natural to take advan- 
tage of the perfect correlations in graph states for the 
constructions of GHZ paradoxes. In this Letter we shall 
identify those graphs, called GHZ graphs, whose corre- 
sponding graph states lead to genuine multipartite mul- 
tilevel GHZ paradoxes. Furthermore we derive a Bell in- 
equality for multipartite and multilevel systems as well as 
a state-independent KS inequality for every GHZ graph. 

As a graph state for qubits is related to a simple graph, 
a nonbinary graph state 27-29] is associated with a 
weighted graph. Let Z d = {0, 1, . . . , d— 1} denote the ring 
with addition modulo d. A Z^-weighted graph G = (V, T) 
is composed of a set V of n vertices and a set of weighted 
edges specified by the adjacency matrix T, a symmetric 
n x n matrix with zero diagonal entries and the matrix 
element T uv G Zd denoting the weight of the edge con- 
necting the vertices u and v. A graph is connected if for 
any pair of vertices u, v there exists a finite number of 
vertices {vi\f =0 such that riiLo 1 ^ViV i+1 ^ with u = vq 
and v = vk ■ 

We denote by D v the degree of vertex v £ V which is 
the sum of the weights of all the edges connecting to v 
and by W the total weight of G which is the sum of the 
weights of all the edges. Explicitly, we have 

d v = r„ (v G V), w = x - T uv (i) 

u£V u,v£V 

A GHZ graph is a connected Z^-weighted graph satis- 
fying (i) the degree of each vertex is divisible by d, i.e., 
D v = mod d, while (ii) the total weight is NOT divisi- 
ble by d; i.e., W ^ mod d. From these two conditions 
it follows immediately that the GHZ graph does not exist 
in odd dimensions and uj w — — 1, where u> = e l ~a~. In 
fact, from the first condition, there is an integer t v such 
that D v = dt v for each v G V, and from the fact that the 
total weight W = dt/2 with t = J2 V <£V^ V 1S an integer, 
since T is symmetric, it follows that if d is odd then t 
must be even and thus W is divisible by d. Furthermore, 
in even dimensions, the total weight W is not divisible 
by d if and only if t is odd and thus u> w — (—1)* = — 1. 
In what follows we shall always assume d to be even. A 
GHZ graph is called " primary" if for each vertex a € V 
there exists a pair of vertices b, c such that r o f, and r oc 
are coprime and "weakly primary" if there exist three 
vertices a,b,c£ V, such that T ab is coprime with T ac . 

In the case of d = 2 a GHZ graph has an odd number 
of edges and every vertex has an even number of neigh- 




FIG. 1: Examples of GHZ graphs. Unlabeled thin black or 
red edges have weight 1 or d — 1, respectively. All possible 
GHZ graphs on 3 and 4 vertices are shown in (I) and (II), 
where a 1 = |+a, b' = |+6, and c' = f+c with a+b+c = d/2. 
A GHZ graph on 5 vertices is shown in (III) where the thick 
red edges have weight d/2 — 1. In (IV) and (V) two primary 
GHZ graphs on 2k + 4 and 2k + 5 vertices (k > 1) are shown. 

bors. All GHZ graphs for d — 2 are primary. For exam- 
ple, a loop graph with an odd number of vertices and a 
complete graph with 4j + 3 (j > 0) vertices are possi- 
ble GHZ graphs. There is only a single GHZ graph on 3 
vertices as shown in Fig. 1(1) and it is clear that it is not 
weakly primary if d > 2. In the case of n — 4 all possible 
GHZ graphs are shown in Fig. 1 (II) with weights satisfy- 
ing a+b+c = d/2. lfd = 4k then d/2±l = 2fc±l are co- 
prime and thus, by choosing, e.g., a = 1, c = 1, we obtain 
a primary GHZ graph. If d = 4k + 2 then there always 
exists a vertex with all edges having even weights, since 
d/2 is odd, so that only a weakly primary GHZ graph 
exists in this case. Examples of primary GHZ graphs 
for arbitrary n > 5 and even dimensions are shown in 
Fig. 1 (III- V). The primary GHZ graph on 5 vertices as 
shown in Fig. 1 (III) can be generalized to any odd num- 
ber of vertices. 

Consider a system of n particles each of which has d 
energy levels, a audit for short, and label them with V. 
Let {|s)„|s € Zd} be the computational basis for qudit 
v G V and {|s)|s G Z^} is a basis for n qudits where Z^ 
is the set of all n-dimensional vectors s = (si, s%, . . . , s n ) 
with components s v G Zd for all v G V. To any weighted 
graph G = (V,T) on |V| = n vertices we can associate 
with a qudit graph state 

l r > = i E ^ S ' r - S |s>, (2) 

« 2 v 
s ez v d 

which is also the unique joint +1 eigenstate of n com- 
muting vertex stabilizers 

g v =X V [] Z£», (3) 

u£V 
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i.e, g v \T) = \T) for all v £ V. Here we have 
introduced the generalized bit shift operator X v — 
J2 S £i, d \( s + 1) m °d d) (s\ v and phase shift operators 
Z v = J2sei.d a;S l s )( s li> f° r eacn qudit v £ V. It is easy to 
check that = Zf, = I and Z V X V = ljX v Z v . Our main 
result reads as follows: 

Theorem For each (weakly) primary GHZ graph G = 
(V,r) on |V| = n vertices, with weights taken values in 
lid, the graph state |T) provides a (weakly) genuine n- 
partite d-level GHZ paradox. 

Before embarking on the proof we should clarify what 
we mean by genuine n-partite and d-level and give an ex- 
ample. According to [IB] a GHZ paradox, formulated via 
a set of commuting observables, is said to be genuinely 
n-partite if one cannot reduce the number of parties and 
still have a Mermin-GHZ paradox. A GHZ paradox is 
(weakly) genuine d-level if one cannot reduce the dimen- 
sionality of the Hilbert space of (all) any one of the par- 
ties to less than d and still have a paradox. 

As an example let us consider the GHZ graph as shown 
in Fig. 1 (II) in the case of n = 4 and the following 5 
commuting observables that stabilize the corresponding 
graph state 

X Z a ' Z b Z c +1 
Z a ' X Z c ' Z b ' +1 
Z b Z c ' X Z a +1 (4) 

Z c Z b' Z a X +1 

A f A f A f A f -1 

which provide us a GHZ paradox. Measurement of the 
product of the operators in each row gives a certainty 
result 1 or —1 as listed in the right column of Eq.Q by 
quantum mechanics. With the analogue to EPR's argu- 
ment, the result m% or m% of measuring the correspond- 
ing d-outcomes measurements X v or Z v on the vth qudit 
can be predicted in advance with certainty with the help 
of the results of spacelike separated measurements of X 
or Z on the other three qudits and are therefore elements 
of reality. Because the algebraic relations are preserved, 
we have m%' z = ui k with u) — e l27T / d for some k £ Id and 

(™?)(m!)°'(mf)Vl) C = 1 

Kr'K)K) c 'K) k ' = i 

K) 6 (mf) c '(mf)(mfr = 1 

(rofr'M-^mirVS)- 1 = -1. (5) 

The contradiction lies in the fact that all five equations in 
Eq. ([5]) cannot hold simultaneously. In the case of d = 4 if 
we choose a = b = 1 and d = 2 with a' = b' = 3 and c = 
then the GHZ graph is primary and the corresponding 
GHZ paradox is genuine 4-partite and 4-level. In the case 
of d = 6 we can choose a = b = c = 1 and a' = b' = d = 4 
such that for the second qudit there exists a projection to 



a qutrit by identification Z 2 for d — 6 with Z for d = 3. 
Thus it provides an example of weakly genuine 6-level 
GHZ paradox that can be regarded as GHZ paradox on 
a hybrid system of three 6-level system plus a qutrit. 

Proof.— Let G = (V,T) be a GHZ graph; i.e., the de- 
gree of each vertex D a is divisible by d and the total 
weight W satisfies lu w = — 1. For each qudit v € V we 
measure two unitary observables X v and Z v with out- 
comes assigned to values mJ;,TO 2 e £ I'd}: respec- 
tively. First of all these values are elements of reality 
because of the perfect correlations g v \T) = \T) (v e V). 
In any local, or noncontextual, hidden variable models 
these values are independent of which observables might 
be measured by other observers. Furthermore, they must 
satisfy the same algebraic rules, e.g., the product rule, as 
their quantum counterparts do. For example from the 
definition of the vertex stabilizer g v it follows 

M v := ml J] (m 2 /"" = 1 (6) 

for each v G V. On the other hand from the constraint 
Ay|r) = — because of the identity 

Y[g a =«, w Xvl[Z?>=-X v , (7) 

a£V aeV 

it follows that rii,ev TO S ~ ~ 1 which is impossible be- 
cause n„ey M v = n«ey m in m which the fact that D v 
is divisible by d has been used. 

By definition a GHZ graph is a connected graph and 
thus for each partition of n observers into two groups 
some of n + 1 unitary observables will not be commuting 
when restricting to either one of two groups. Therefore 
the GHZ paradox for |T) is a genuine n partite. Further- 
more, if the GHZ graph is primary then each vertex is 
attached to at least one pair of edges of coprime weights. 
If there were a projection to lower dimensions for a qudit, 
some eigenstates of X a and those of Z^ ab and Z v a ac are 
orthogonal. This is impossible because first there always 
exist p,q £ Zd such that pT a b + qT ac = 1 mod d and, 
second, X v and Z v are two complementary observables 
whose eigenstates cannot have a zero overlap, which is 
the case if the dimensionality can be reduced. (t 

Some remarks are in order. First, for we have con- 
structed genuine n-partite and ci-level GHZ paradox with 
n > 5 can be even. Second, any state that is related with 
GHZ graph states via local unitary transformations ex- 
hibits also GHZ nonlocality. Third, for a graph that is 
not GHZ graph it is also possible to construct a GHZ 
paradox for the graph state if the underlying graph con- 
tains a GHZ subgraph. A subgraph H = (V',V) of a 
weighted graph G = (V,T) is also a Z^-weighted graph 
with a vertex set given by V C V and edges specified 
by T' ab — T ab if a, b £ V. If furthermore the subgraph 
is a GHZ graph we shall refer to it as a GHZ subgraph 
of G. Suppose that the graph G contains a GHZ graph 
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H = (V, T') with \V'\ = m < n, then the m + 1 observ- 
ables g u with u £ V and Iluev 9u yield a GHZ paradox 
for the graph state |T) . It is clear that it is only a genuine 
TO-partite GHZ paradox if the GHZ subgraph is primary. 
For example, the 4-qubit GHZ state is equivalent to the 
graph state corresponding to the complete graph on 4 
vertices, which contains a loop of length 3 as a GHZ 
subgraph. In fact the original GHZ proof [4] revealed a 
3-partite GHZ nonlocality using this GHZ subgraph. 

As the first application we shall derive a Bell inequality 
with two measurement settings for each observer with 
the help of the GHZ paradox derived from a GHZ graph. 
Consider two d-outcome measurements A v and B v for 
each observer v G V and assign values in £ Z<j} to 
them (Bell-KS value assignment). For each GHZ graph 
G = (V, r) we introduce a Bell operator as 



Ba = 




B 



AT,, 



(8) 



v£V 



uev 



v£V 



Taking into account the identity X)fc=o = d$l,o f° r 
arbitrary I £ and denoting Ay — Y[ v ev an d N v = 
A v Yiuev Bu" v f° r each v £ V, where Si.j is the standard 
Kroneckcr delta symbol, we have 



Bo = S. 



l, A x 



-£(*. 



N„ 



S-i,n v ) < n-1. (9) 



The inequality holds in any local realistic theory because 
if there are n positive terms then there is necessarily a 
negative term in Bq- If N v = 1 for every v £ V then it 
holds Ay = 1 which contributes a negative term; if N v — 
1 for all v £ V — {v } and Ay = — 1 then it necessarily 
holds N Vo = — 1 because Ay = Y\ veV N v . Furthermore 
it is easy to see that Bq < n + 1, which is attained by the 
graph state |T) with (r|Sc|r) = n+1 in which A v and B v 
are chosen to be X v and Z v , respectively, for each v £ V. 
In this case the quantum to classical ratio (n + l)/(n— 1) 
is a constant independent of the dimension, comparing 
to that of [3U]. 

Every GHZ paradox leads also to a proof of KS the- 
orem. And any proof of KS theorem can be converted 
to an experimentally testable inequality, called as KS in- 
equality, in the manner of Cabello [8]. As the second 
application we consider the following KS inequality 

~ /xi n Xv + e aix v n z i uv + h - c ) 

\ v£V vGV uGV I c 

-2(4n 3v + H.c.J < C n+X4 (10) 

where, with A = 2 (n+i) arLC ^ ^ = 27r/d, we have denoted 
= (A- LAJ)cos[A]6»+(l + [AJ - A)cosLAJ0. (11) 



First, each term, e.g.,(^X v Y\ v€V g v ^j , is the abbreviated 

form of the classical correlation of n+ 1 observables, e.g., 

([x v gig2-..gi^J ■ Second, the upper bound can be easily 

inferred from the Lemma proved below. Third, we have 
C„+i i( 2 < n + 2 while the quantum mechanical value of 
the left-hand side of Eq.(10) equals to n + 2 identically 



and therefore violates the above KS inequality in a state- 
independent fashion. 

In summary, first of all we have identified a special 
kind of graphs, called GHZ graphs, whose correspond- 
ing graph states give rise to GHZ paradoxes. Except 
for the case n = 4 with d = 4k + 2 for which only a 
weakly genuine GHZ paradox is found we have derived 
genuine n-partite and eHevel GHZ paradoxes from qudit 
graph states corresponding to GHZ graphs with n > 4 
and even d being arbitrary. Second, as applications for 
each GHZ graph we derive a Bell inequality with two 
d-outcome observables for each observer whose maximal 
violation is attained by the corresponding graph state as 
well as a state-independent KS inequality that is satis- 
fied by any noncontextual hidden variable models. This 
would be helpful to the analysis of multipartite contex- 
tuality or multipartite nonlocality. It should be noted 
that GHZ paradoxes may exist for those states that are 
equivalent to the graph states under local Clifford (LC) 
transformations. However the conditions under which 
both two GHZ paradoxes arising from two LC equiva- 
lent states are genuine n-partite seem to lie out of the 
reach of current Letter. Besides, the examples we are 
analyzing here involve only some special classes of graph 
states, so figuring out other classes of graph states which 
are consistent with our theorem are still meaningful as 
for fixed parties n, different graphs may have different 
robustness against decoherence, which may help to de- 
sign new quantum protocols for reducing communication 
complexity Ironically, a genuine 4-partite GHZ paradox 
is still missing for the original 4-qubit GHZ state. 
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Lemma Let = {9, 29, . . . , d9} with 



A 



and x 



2(n+l) 

variables and x 



max < fix) :- 



(x 1 ,X 2 , ...,X n ) 

Xi. We have 



2n/d and 
be n real 



i=l 



cosxi - cos (x s ) \ = C„ ;d . (12) 



Specially, if n > d/2, i.e., A < 1, then C n d = n + 1 — 
dsin 2 ^ and if d = 2{n + 1)1 for some I, i.e., A is an 
integer, then C nd = (n + 1) cos 



Proof. The maximum of ,f(x) over W 1 is the largest 
value on all local extremal points satisfying = 
sa\x s — sinxi = (Vi). Let x\ = x then cither Xi = 
a,i = 2uiTT + x or x.i = bi = (2ui + l)ir — x for all in- 
tegers Ui with i > 2 since sinx^ = sinx. Denote by m 
the number of Xi's being equal to <Xj and k = n — m 
the number of x^s being equal to frj's among {xi}™ =l . 
Then from sin a; = sin(fc7r + (m — k)x) it follows either 
a) x = x a with 2ln + x a = kir + (m — k)x a or b) x = X\> 
with (21 + 1)tt — Xb = kir + (m — k)xb for all the in- 
tegers I > 0. At these extremal points we have either 
f(x a ) = (to — fe — 1) cos x a or f(xb) = (to — k + 1) cos x&. 
If k > 1 then m — k = n — 2k < n — 2 and thus 
f{xa,b) < n — 1. If k = then f(xi) = (n + l)cosa;j, 
where x\ = (xi, xi, ... , xi) with xi — (21 + l)ir/(n + 1). 
Since f(xo) > n—l and f(xo) > f(x*i) the extremal point 
xq leads to the largest value of /. 

If A is an integer then xq € and thus the global max- 
imum f(xo) is attainable and in this case f(xo) — C n d- 
If A is not an integer then the maximal value f(xo) is 
not attainable by any vector in U^- However its maxi- 
mum must be attained at those vectors near one of those 
extremal points that have the floors or ceilings of the 
components of the extremal points. We consider at first 
those vectors in 11% near xq that have a number m of 
x + = [A] 9 and a number n — m of x_ = [AJ 9 as com- 
ponents with n > m > 0. On these vectors f(x) assumes 
values F m — m cos x + + (n — m) cos x_ — cos(m9 + 
Let A TO = (F m+ i — F m )/(2sin |) and we have 



. 2nx_ + (2m + 1)9 . x_ + x + 
sin sin 



(13) 



Since A is not an integer we have < (x_ +x + )/2 < ir/2. 
In the case of d/2 > n we have nx- + m9 < 2tt + x_ and 
then A m > if m < S := d/2 - (n + 1) |AJ and A m < if 
m > 5. As a result maxF m = Fg = C n d- If n > d/2 then 
x_ = and in this case A m > if m' u < m < m u with 
m' u = ud or m u = (u + l/2)d — 1 and A m < otherwise 
for u > being integer. Thus the maximum value of F m 
must be taken on m' u or m u and obviously F mo > F mu 
and F ma > F m '^ for all u. As a result we have maxF m = 
F mo = Cnd- Since C nd > (n + l)cos[A]0 > f(x t ) for all 
I > 1 we see that C n d is the global maximum of /. ft 



